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Abstract 

We study a simple model of DNA divalent cationic surfactant complexation. We 
find that the combination of electrostatic and hydrophobic effects leads to a coop- 
erative phenomenon in which as the amphiphile is added to the solution containing 
DNA, a large fraction of the DNA's charge is neutralized by the condensed diva- 
lent cationic surfactants, forming the surfoplex. This binding transition occurs for 
concentrations that are lower for divalent than for monovalent surfactants. Since 
the electrostatic strength is larger in the first case and the amount of surfactant 
lower, we suggest that multivalent amphiphilic molecules would be more efficient 
than monovalent for transfection. 
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1 Introduction 



Gene therapy represents a promising way for the treatment of both genetic 
and acquired diseases. The basic idea is to replace the sick gene with a healthy 
one or in some cases to add a new gene to get a resulting synthesis of a thera- 
peutic protein. The process of getting the new gene to its target involves the 
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crossing of several barriers, such as the cell membrane and the nuclear mem- 
brane. Since both the DNA and the cell membranes are negatively charged, 
the naked polynucleotides are electrostatically prevented from entering into 
the cell. Viral vectors such as retroviruses and adenoviruses are very efficient 
and able to target a wide range of cells [1] [2] . The gene transfection is accom- 
plished by the use of a virus in which the native DNA has been replaced by 
the required DNA. Since the main purpose of the virus is to replicate itself, 
the new gene is successfully transfected to the cell nucleus by endocytosis or 
membrane fusion. This process suffers from the drawback that in cases where 
repeated treatment is needed, this might cause the immune system to react 
negatively due to the viral origin of the vector. 

As a solution for this problem nonviral vectors have been developed. One of the 
approaches, pioneered by Feigner and Ringold [3] [4] relies on the association 
between anionic nucleic acid and cationic lipid lipossomes. The process of 
association neutralizes the excess negative charge of the DNA and the DNA- 
lipossome penetrate into the cell by endocytosis. The efficiency of this method 
is low and they are toxic to the cell at the concentrated used. 

Among the parameters needed for achieving efficient transfection is the re- 
quirement that the complex should have small size similar to a virus [5] . This 
factor indeed limits the efficiency of the DNA-lipossomes complexes. In con- 
trast, cationic surfactant have been shown to condense to DNA into discrete 
particles containing a single nucleic acid molecule [6] [7] that can neutralize 
and revert the charge of the DNA [8]-[ll] allowing the surfoplex to approach 
the membrane. Monovalent DNA-surfactant complexes exhibit a discrete first- 
order phase transition between elongated coil and collapsed globule [12] for 
a concentration of amphiphilic molecules well below the miscclar concentra- 
tion. Despite this unique feature, monovalent detergent are poorly efficient in 
vitro in gene transfer [13] [14]. When the surfoplex approaches the membrane, 
the interaction between the surfactant molecules and the phospholipid bilayer 
overcomes the electrostatic attraction between the anionic nuclei acid and the 
cationic detergent. The surfactant incorporate into the phopholipid membrane 
resulting in the unfolding of the DNA that stays outside the cell [15]. 

Therefore, in order to have a stable surfoplex, one needs to enhance the elec- 
trostatic interaction. In this paper we present a model of DNA-amphiphilic 
solution where the surfactant is divalent. We find that in equilibrium, solu- 
tion consists of complexes composed of DNA and associated counterions and 
amphiphiles. Even for an small amount of amphiphiles in the solution, the co- 
operative binding is found. Due to the high valence, the interaction between 
the nuclei acid and the detergent is expect to overcome the interaction be- 
tween the surfactant molecules and the phospholipids in the membrane and 
stay associated to the DNA allowing it to transpass the membrane. The com- 
plex formed is more stable than the monovalent structure, forming a collapsed 
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Fig. 1. The model 
globule that we expect will transect with high efficiency [16]. 

2 The Model 

Our model, illustrated in Fig. 1, consists of a solution of DNA strands of 
length L and diameter a^, divalent surfactant and monovalent salt. In aque- 
ous solution, the polyions become ionized resulting in a negative charge —Zq 
distributed uniformly distributed with separation h = L/Z. The solvent, wa- 
ter, is modeled as a continuous medium of dielectric constant D. The ions of 
the salt are completely dissociated, forming an equal number of positive and 
negative ions. Similarly, the surfactants are assumed to be fully dissociated 
producing negative monovalent colons and polymeric chains with a divalent 
cationic head group. For simplicity, all the counterions and colons are treated 
as identical, independent of the molecules from which they were derived. The 
electrolytes are depicted as hard spheres of diameter Oc and charge ±g and 
the surfactant is modeled as a polymer of s monomers each one considered 
as a rigid sphere of diameter with the head monomer carrying a charge of 



The interaction between the hydrophobic tails is short ranged and character- 
ized by the parameter x- The density of DNA strands is pp = Np/V, the den- 
sity of monovalent salt is Ps — N^/V and the density of divalent amphiphiles 



The strong electrostatic interactions between the polyions, the counterions, 
salt cations and surfactant heads leads to the formation of complexes, which 
in thermodynamic equilibrium will be made up of one polyion, ric monovalent 
counterions and divalent surfactant. We do not consider the effects of of 
polydispersity in the size of the complexes, since it does not affect the final 



+2q. 



is Pa = NJV. 
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result. Due to the association and to the charge conservation, there are only 
two free quantities and so, 

Pc^{Z - nc)Pp + Ps Pa+ ^ Pa- riaPp P- ^ Ps + 2pa (1) 

where pc is the density of free monovalent counterions, pa+ is the density of 
free amphiphiles and p_ is the density of negative ion. 

The objective of this theory is to determine the number of counterions Uc 
and surfactants associated to each DNA strand. For this, we construct 
the Helmholtz free enery of the system and minimize it. The details about 
the model can be found elsewhere [8]- [11]. We give here the main steps. The 
relevant contributions for the Helmoltz free energy are two, the electrostatic 
and the entropic namely: 

F^Fel + Fent (2) 



In Fci three types of interactions can be found: between the free ions and free 
surfactants ,Fis, between the complex, free ions and free surfactants ,Fpis, and 
between the complexes Fpp. 

Fei F-ig -\- Fpig -\- Fpp (3) 



With the aid of the theory of Debye-Hiickel-Bjerrum (DHBj) it is possible to 
find the electrostatic interaction between the complexes, ions and surfactant 
given by [8]- [11]: 



with (5 = l/ksT and 

- / 

^ ' 



where k in (/ta)^ = Anpl/T* is the inverse of the Debye screening length, 
pI = pia^ = PcO? + 4pa+a^ + P-i^ is the reduced density and T* = DkBTa/(f 
is the reduced temperature. Furthermore, Zc — Z — Uc — 2na is the valence 
of each complex and a — {ac + dp)/ 2 is the effective radius of the exclusion 
cylinder around each complex. 

In the framework of the Debye-Hiickel theory, the interaction between the free 
ions is given by: 
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ln(l + uttc) — uttc + 



[nar 



(6) 



The electrostatic free energy interaction between two complexes for large sep- 
arations is screened. The short-range of this interaction allows one to use a 
mean-field approximation resulting in: 

^ 27ra3Z,V^exp(-2/ta) 



The calculation of the entropic contribution can be obtained with the aid of 
Flory [17] theory of mixing. The free energy is a sum of ideal free energies of 
various species, namely 



ent 



E 



Ps - Ps In 




(8) 



where s represents the different species and is the internal partition of the 
species s. In the case of the particles without structure, the internal partition 
function C,_ — C,^ — C,^_^ — 1. The volume fraction 0^ of the different species 
are: 




Here we introduce the fractions counterions, rric = Ua/Z, and surfactant, — 
Uc/Z, associated to each DNA strand. Accounting for all these terms, the 
entropic contribution becomes 



+ Pcln 0c - Pc + Pa+ In ^ - 



a 



+ p_ln0_-p_ (10) 

where the entropy for the free surfactants arises from the Flory theory for poly- 
mers [17] which is also the basis of the entropic contribution for the complex 
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[8]. The internal partition function of the complex. (,.1- can be calculated by 
modeling the DNA by an one dimensional lattice with Z sites. If the number 
of associated ions to each site can be only zero or one, this problem becomes 
equivalent to finding the free energy of an one dimensional array with the 
three different states. 

Due to the presence of two different valences, the exact solution of this model 
is not trivial, so we employ the mean-field Gibbs-Bogoliubov-Feymman in- 
equality. The resulting partition function is given by 



InCc [rric, ma] = 



y 2 

^ -1 

Z2 



+ /3X(^-I)m^ 



-f ZrUc In mc + Znia In nia 

+ Z{l-mc- ma) In (1 - mc - m^) 



(11) 



where ^ = I3q^ / Da is the Manning parameter, K — Z[ip{Z) — '0(1)] — Z + 1, 
and ip{n) is the digamma function. 



The equilibrium configuration of the system is found by the minimization of 
the Helmoltz free energy, leading to two equations, namely 



Smc = 



dmc 
dF 

dma 



(5ma = 0. (12) 



Solving this system of two equations, it is possible to obtain the values of mc 



and ma- 



3 Results and Conclusions 



We define a "surfoplex" to be a complex in which almost all of the DNA's 
phosphate groups are neutralized by the associated surfactant molecules. As 
mentioned earlier, we are interested in the minimum amount of cationic sur- 
factant needed to transform naked DNA into surfoplexes. To this effect, we 
study the dependence of the number of condensed surfactant molecules on the 
bulk concentration of surfactant Pa- In order to evaluate the relevance of hy- 
drophobic interactions between the amphiphiles, the hydrophobic parameter 
was varied from to I3x = —6. The effect of addition of high concentrations 
of salt to the system was analyzed by varying the amount of salt added to the 
system p^. 
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Fig. 2. Effective binding of amphiphiles 
(a)ma and countcrions (b) rric as a func- 
tion of tlic surfactant concentration pa- 
The density of DNA and salt are respec- 
tively 2 X IQ-^M and lO'^M. The hy- 
drophobicity are Px = —3.5 (solid line) 
and Px = ~6 (dashed line). 

Fig. 2 and Fig.3 illustrate the surfactant and the counterion binding isotherm, 
nia and rric respectively, as a function of total amphiphilic concentration for 
two different hydrophobicity parameter f3x = —3.5, —6 and salt concentrations 
Ps = 10^'^ Af, lO^^M. At very low densities of divalent surfactant, the conden- 
sation is dominated by the monovalent ions, since the divalent amphiphilic 
molecules gain more entropy, and thus lower the total free energy, by staying 
free. However, as the concentration of divalent surfactant increases, the gain 
in electrostatic energy due to condensation and to the hydrophobicity of the 
amphiphilic molecules wins over the entropy. 

We note that unlike the association with the ionic monovalent surfactants, 
which exhibits a large degree of cooperativity characterized by the sharp rise 
in the surfactant binding fraction, the replacement of the condensed mono- 
valent counterions by the divalent surfactant proceeds more smoothly. This 
result could have been anticipated a priori. After the first amphiphile is asso- 
ciated, the condensation of additional molecules is energetically favored since 
the buildup of the hydrocarbon density in the vicinity of a polyion helps 
to exclude water and, thus, reduces the unfavorable hydrophobic energy of 
the alkyl tails. This effect competes with the electrostatic repulsion between 
the like-charged counterions that tends to inhibit it any further association. 
The repulsion is stronger for divalent ions than for monovalents and this ex- 
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Fig. 3. The same as Fig.2 for ps = lO^^^M. 

plain why the cooperative effect is stronger in the later. In the case of the 
divalent amphiphilic molecules the surfoplex is formed for a density of sur- 
factant around 0.0005M as illustrated in Fig. 4 what is much lower than the 
one observed in the monovalent case [8] -[9]. For comparison with the case of 
monovalent surfactant illustrated in refs. [8] for producing Fig. 4 we employ 
the hydrophobicity factor /3x = —4. 

Fig. 5 and Fig 6 illustrates the effective charge of the complex. Due to the 

simplicity of our model that does not allow for the association of more than 
one surfactant molecule to each charged group along the DNA, there is no 
charge inversion. However, the decrease in the charge produced by the divalent 
surfactant is larger than the one observed when the amphiphilic molecules are 
monovalent for a similar model [8]- [9]. 

In resume, we have presented a simple theory of DNA, for monovalent salt 
and divalent surfactant solutions. Our results should be of direct interest 
to researchers working on the design of improved gene delivery systems. In 
particular, we find that addition of cationic divalent surfactants leads to a 
cooperative binding. This binding happens far below the critical micell con- 
centration and far below the concentration for this transition to happen in 
the presence of monovalent surfactant. Until now experimental attempts in 
employing surfactants as nonviral agent for transfection have been limited to 
the use of monovalent surfactants. The surfoplex produced with monovalent 
amphiphilic exhibit poor eficiency. Close to the membrane, most of the surfac- 
tants disassociate from the DNA and the transfection does not occurs [13] [14]. 
We propose the use of divalent cationic surfactants in the formation of the 
complex. Besides being more strongly connected to the charged groups of the 
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Fig. 4. Effective binding of divalent sur- 
factant nia = na/Z and countcrions 
n^c = nc/Z as a function of the surfac- 
tant concentration pa for varius salt con- 
centrations: (a) 5mM, (b) 18mM and (c) 
40mM. The density of DNA is 2 x W'^M 
and the hydrophobicity is f3x = — 4. 

DNA what will give more stability to the complex during the transfection, the 
amount of surfactant required is lower. Since the surfactant are toxic to the 
organism, this should reduce the risk of unnecessary medical complications. 
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